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We present results of large-scale numerical simulations of the evolution of neutrino and antineu- 
trino flavors in the region above the late-time post-supernova-explosion proto-neutron star. Our 
calculations are the first to allow explicit flavor evolution histories on different neutrino trajectories 
and to self-consistently couple flavor development on these trajectories through forward scattering- 
induced quantum coupling. Employing the atmospheric-scale neutrino mass-squared difference 
(|(5m^| ~ 3 X 10"'^ eV'^) and values of 0i3 allowed by current bounds, we flnd transformation of 
neutrino and antineutrino flavors over broad ranges of energy and luminosity in roughly the "bi- 
polar" collective mode. We find that this large-scale fiavor conversion, largely driven by the flavor 
off-diagonal neutrino-neutrino forward scattering potential, sets in much closer to the proto-neutron 
star than simple estimates based on flavor-diagonal potentials and Mikheyev-Smirnov-Wolfenstein 
evolution would indicate. In turn, this suggests that models of r-process nucleosynthesis sited in the 
neutrino-driven wind could be affected substantially by active-active neutrino flavor mixing, even 
with the small measured neutrino mass-squared differences. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

In this paper we employ large-scale computational 
techniques to tackle the vexing problem of neutrino flavor 
transformation in the core collapse supernova environ- 
ment. Neutrinos and the weak interaction play pivotal 
roles in the core collapse/explosion phenomenon. The 
Chandrasekhar mass core of iron-peak material left at 
the end of the hydrostatic evolution of a massive star 
goes dynamically unstable and collapses in ~ 1 s to a 
proto-neutron star configuration at nuclear density. The 
amount of gravitational energy promptly converted into 
trapped seas of neutrinos is ~ 1% (~ 10^^ erg) of the core 
mass. Within a few seconds after bounce ~ 10% (~ 10^'^ 
erg) of the core mass (the gravitational binding energy) 
will be emitted as neutrinos. 

Nearly all of this gravitational energy is converted into 
seas of Ve, '^e , Vfi, i^/j , '^r and v-r neutrinos in rough energy 
equipartition. Though these neutrinos diffuse with short 
mean free paths in the proto-neutron star, they decouple 
near the stellar surface where the matter density falls off 
steeply, the so-called neutrino sphere. Neutrinos propa- 
gate nearly coherently above this point, though neutrino- 
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matter interactions, especially the charged current cap- 
ture reactions v^. + n —^p + e~ and D^+p n + , can 
deposit energy and set the local ncutron-to-proton ratio, 
n/p. 

For this reason, and because the fluxes and energy 
spectra may be different for Ve, i^e and fpZ/pZ^,-P^, the fla- 
vor content of the neutrino field above the proto-neutron 
star and its evolution in time and space can be impor- 
tant [H, 0, Q • This can be true both for the supernova 
shock reheating epoch (where the time post core-bounce 
is tpB ^ 0.5 s) and in the later hot bubble, neutrino- 
driven wind epoch (ipB 3 s) . In this paper we concen- 
trate on the latter epoch. 

Following the development of neutrino and antineu- 
trino flavors in the coherent regime above the proto- 
neutron star surface is challenging. The potential gov- 
erning the effective neutrino mass differences in this en- 
vironment will have contributions from charged current 
neutrino-electron forward scattering and neutral-current 
neutrino-neutrino forward scattering. The former contri- 
bution 4] is diagonal in the flavor basis, while the latter 
neutrino-neutrino potential has both flavor-diagonal ^5] 
and flavor off-diagonal [1, 0] components. The neutrino- 
neutrino forward scattering potential renders the neu- 
trino flavor evolution problem nonlinear in the sense that 
the potential which governs neutrino flavor transforma- 
tion is itself dependent on the flavor evolution histories 
of the neutrinos. 

Furthermore, neutrinos propagating on intersecting 
world lines can have their flavor evolution subsequently 
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quantum mechanically coupled by forward scattering (see 
Fig. 8 in Ref. (3|] and the text beneath it). We some- 
times will refer to this coupling as "entanglement" . By 
this terminology we do not mean quantum entanglement 
of momentum states, a phenomenon which has been ar- 
gued to be unimportant in the supernova environment Q . 
In any case, neutrino trajectories coming off the proto- 
neutron star surface at different angles in general will 
have different flavor evolution histories which must be 
self-consistently calculated. 

Another issue revolves around the efficacy of a mean 
field Schrodinger-like or Boltzmann kinetic equation ap- 
proach to the evolution of neutrino flavors [1, H, [^, [l^ . 
In this paper we will assume that higher order correla- 
tions in neutrino-neutrino scattering are unimportant in 
the coherent regime above the proto-neutron star. 

Previous attempts to model neutrino flavor evolution 
in the coherent regime P, H, H, H, [ll|, [H, [l^ have made 
the approximation that all neutrinos evolve in flavor 
space the way a radially-propagating neutrino does. This 
we will term the "single-angle" approximation. Since the 
neutrino-neutrino forward scattering potential is inter- 
section angle dependent, this is not always a good ap- 
proximation, especially for regions close to the proto- 
neutron star. 

However, these previous studies done with the single- 
angle approximation have found that it is possible to have 
large-scale collective behavior in neutrino flavor evolu- 
tion, where all, or some significant subset of, neutrinos 
experience similar time/space flavor evolution histories. 
They also have shown that neutrino flavor transforma- 
tion can differ significantly from the Mikheyev-Smirnov- 
Wolfenstein (MSW) 0, Q [ll paradi gm. Recent work 
[TgI has shown that the expected neutrino fluxes in both 
the shock reheating and hot bubble epochs could provide 
the "necessary" conditions for large-scale simultaneous 
collective neutrino and antineutrino flavor transforma- 
tion over broad ranges of neutrino energy. Whether these 
expected neutrino fluxes are actually "sufhcient" to ob- 
tain these collective modes has remained an open ques- 
tion, to be answered with an appropriately sophisticated 
numerical simulation. Likewise, the range of possible col- 
lective neutrino behavior be it the "synchronized" 
mode ,18] or the "bi-polar" mode [l3], niay depend sen- 
sitively on the neutrino flux conditions and on the geom- 
etry. 

It should be noted that many previous numerical stud- 
ies employing the single-angle approximation have also 
used unphysically large values of neutrino mass-squared 
difference. This is because with straight MSW, and with- 
out taking account of the neutrino-neutrino scattering- 
induced flavor off-diagonal potential, it requires \Sm'^\ > 
1 eV^ to have significant neutrino flavor transformation 
deep enough in the supernova envelo pe to affect shock 
reheating or the r-process [U, H, H, [H, HO]- 

However, recent observations/experiments (see, e.g., 
Ref. [2l[ for a review) have revealed much about the fun- 
damental flavor mixing parameters of the three known 



"active" neutrinos. (In this paper, we will ignore the ef- 
fects of speculative additional "sterile" neutrino states.) 
We know the two mass-squared differences, the atmo- 
spheric scale, <5mat,n — 3 x 10~^eV^, and the solar scale, 
SniQ ~ 8 X 10~^ eV^. We as yet do not know the neutrino 
mass hierarchy related to the atmospheric mixing and 
we do not know the absolute neutrino mass eigenvalues. 
Of the four mixing parameters in the unitary transfor- 
mation between the flavor (weak interaction) eigenstates 
and the mass (energy) eigenstates, we know two of the 
three vacuum mixing angles, O12 and ^23, and we have a 
firm upper limit on ^13, sin^ 26'i3 J; 0.1. We do not know 
the CP-violating phase. 



In this paper we study 2x2 neutrino and antineutrino 
flavor transformation at the Sm^^^ scale, explicitly fol- 
lowing the coupled flavor evolution on neutrino trajecto- 
ries ranging from radially-directed to those tangential to 
the neutron star surface. (In other words, we perform 
"multiangle" calculations with many trajectory/angle 
zones.) Our goal is to study the nonlinear behavior of 
the neutrino field in the coherent regime and to find 
out if large-scale (collective mode) neutrino/ antineutrino 
transformation can occur in the late-time supernova envi- 
ronment. We specialize to late time for two reasons: (1) 
this epoch is when there may be significant differences in 
flux or energy spectrum between Ug, Dg, and the mu and 
tau flavor neutrinos; and (2) this epoch may have a sim- 
pler, more compact, matter density profile near the neu- 
tron star surface. We follow Refs. [22, [2^ and argue that 
2x2 mixing is adequate because the and v-r neutri- 
nos are nearly maximally- mixed in vacuum {623 — 7r/4) 
and these species experience nearly identical interactions 
everywhere in the late-time supernova environment. 



In Sec. we summarize the physical and geometric 
assumptions in our numerical simulations in what we 
call the "neutrino bulb model". In this section we also 
present the basic physics of neutrino flavor transforma- 
tion in the practical formalism used in our numerical sim- 
ulations. We also review the spin analogy for neutrino 
flavors, and estimate the (MSW) resonance locations in 
the hot bubble using both the standard MSW and syn- 
chronization mechanisms. In Sec. IIIII we explain some 
details of our numerical codes and discuss the particular 
numerical difficulties and potential pitfalls in multian- 
gle simulations. We also present the main results of our 
multiangle simulations. The simulations show large-scale 
flavor transformation different from what would be pre- 
dicted if the conventional MSW or synchronization mech- 
anisms apply. In Sec. II VI we identify the flavor transfor- 
mation in our results as being of the bi-polar type [l3|, 
and we analyze this behavior with the help of single-angle 
simulations. In Sec.|V]we give our conclusions. 
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II. BACKGROUND PHYSICS 



A. Neutrino Bulb Model 

At tpB ^ 3 s, the inner core of the progenitor star has 
settled down into a proto-neutron star with a radius of 
about 10 km. In the following ~ 10 s, the nascent neu- 
tron star radiates away its gravitational binding energy 
as outlined above. During this time, neutrinos could de- 
posit energy into the matter above the neutron star and 
create a high-entropy "hot bubble" between the proto- 
neutron star surface and the shock. Inside the hot bub- 
ble, a quasistatic and near adiabatic mass outflow, the 
so-called "neutrino-driven wind" , may be established at 
this epoch as a result of neutrino/ antineutrino heating 
[U [1^. To simplify the numerical calculations of the 
flavor transformations of neutrinos and antineutrinos in- 
side the hot bubble, we approximate the physical and 
geometric conditions of the post-shock supernova by a 
"neutrino bulb model" . This model is characterized by 
the following assumptions: 

1. The neutron star emits neutrinos uniformly and 
isotropically from the surface of a sphere (neutrino 
sphere) of radius R^; [Note that the neutrino flux 
emitted at angle do with respect to the normal di- 
rection at the neutrino sphere comes with a geo- 
metric factor cosi^o- See Eq. ((S]).] 

2. At any point outside the neutrino sphere, the phys- 
ical conditions, such as baryon density nb, temper- 
ature T, etc., depend only on the distance r from 
this point to the center of the neutron star; 

3. Neutrinos are emitted from the neutron star sur- 
face in pure flavor eigenstates and with Fcrmi-Dirac 
type energy spectra. 

The neutrino bulb model, as illustrated in Fig. [1] has 
multifold symmetries. It is clearly spherically symmetric. 
This means that one only need study the physical condi- 
tions at a series of points along one radial direction, which 
we choose to be the z-axis. It is also obvious that the 
neutrino flux seen at any given point on the z-axis has 
a cylindrical symmetry. As a result, different neutrino 
beams possessing the same polar angle with respect to 
the z-axis and with the same initial physical properties 
(flavor, energy, etc.) should be completely equivalent. 
In other words, they will have identical flavor evolution 
histories. One may choose this polar angle to be d, the 
angle between the direction of the beam and the z~axis. 
Alternatively, a beam could be specified by the polar an- 
gle Q giving the emission position of the beam on the 
neutrino sphere (see Fig. [T]). A third option, which we 
have found to be most useful in our numerical calcula- 
tions, is to label the beam by emission angle ^q. This is 
defined to be the angle with respect to the normal direc- 
tion at the point of emission on the neutrino sphere (see 
Fig. [1]) . This emission angle z?o is an intrinsic geometric 
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FIG. 1: The geometric picture of the neutrino bulb model. An 
arbitrary neutrino beam (solid line) is shown emanating from 
a point on the neutrino sphere with polar angle O. This beam 
intersects the «-axis at point P with angle Because neutri- 
nos are emitted from the neutrino sphere of radius 7?^, point 
P sees only neutrinos traveling within the cone delimited by 
the dotted lines. One of the most important geometric char- 
acteristics of a neutrino beam is its emission angle i^o, defined 
with respect to the normal direction at the point of emission 
on the neutrino sphere (t^o = -)- i9). All other geometric 
properties of a neutrino beam may be calculated using radius 
r and i^o- 



property of the beam, and does not vary along the neu- 
trino trajectory. Moreover, because of assumptions[l]and 
[2]in the neutrino bulb model, all the neutrino beams with 
the same emission angle i?o and the same initial physical 
properties must be equivalent. In simulating the flavor 
transformations of neutrinos in the neutrino bulb model, 
it is only necessary to follow a group of neutrinos which 
are uniquely indexed by their initial flavors, energies and 
emission angles. 

At any given radius r, all the geometric properties of 
a neutrino beam may be calculated using r and "d^- For 
example, d and 8 are related through the following 
identity: 



sin d sin G sin i9n 



Rv 



I -In 



where 



and 



I = r cos d, 



Iq = Rv cosi?o- 



(1) 



(2) 



(3) 



Length I — Iq in Eq. ^ is also the total propagation 
distance along the neutrino beam. At a point at radius 
r, the neutrino beams are restricted to be within a cone 
of half-angle 



i^max = arcsin 



Rv 



(4) 



(see Fig. [J). 

One can integrate flux over all neutrino beams (angles) 
and calculate the neutrino number density at radius 
r. In this paper we use the symbol v in the general 
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sense, denoting either a neutrino or an antineutrino. We 
use (va) to denote a neutrino (antineutrino) in flavor 
state a, and Va [va] to denote a neutrino (antineutrino) 
created at the neutrino sphere initially in flavor state a. 
As an example, we shall calculate the differential number 
density dn^^ (q) at radius r: this will have contributions 
from all with energy q which propagate in directions 
within the range between q and q + dq. Here a hatted 
vector n denotes the direction of vector n, and is defined 
as n = n/|n|. The differential number density drtp^(q) 
of can be calculated in a similar way. One finds that 



.(q) = 



jy^ [q] COS doRl d(cos 6)d<i> 



{I - kr 

= ji^ciq) d{cos-d)d(f), 



(5a) 
(5b) 



where the velocity of a neutrino is taken to be the speed 
of light (c = 1), and j^^ (q) is the number fiux of Va with 
energy q emitted in any direction at the neutrino sphere. 
In Eq. (|5ap . i?^ d(cos 0)d$ is the differential area on the 
neutrino sphere which emits neutrinos in the directions 
within the range between q and q+dq, and the factor {I— 
Iq)^'^ accounts for the geometric dilution of the neutrino 
density. In Eq. (|5bp , and are the polar and azimuthal 
angles of q and, in deriving the equation, we have used 
Eq. ([T]) and the identities 



d$ = d<j), 
cosi?oi?,.de = (/-/o)di?. 



(6) 
(7) 



As an added check on Eq. ([5]), note that the total num- 
ber of i>a with energy q passing through the sphere of 
radius r per unit time is 

47rr^ J cos z9 dn^^ (q) 

= 8Tr^r^jy^{q) f cosi?d(cosz?) (8a) 

J COS l9max 

= An^RljyJq). (8b) 

This is indeed equal to the number of i^a with energy q 
emitted per unit time from the neutrino sphere, 

AttRI f 2^>^(g)cosi?od(cos?9o) = 4^'i?'i.^(g). (9) 
Jo 

We note that this flux can also be expressed as 

"1 



A-KRf, / 27rjy^(g)cos7?od(cos?9o) 



where L^^ , {E^^ ) and f^^ (q) are the energy luminosity, 
average energy and normalized energy distribution func- 
tion of Va , respectively. Therefore one has 



ii^c (9) 



(11) 



We take fi,{q) to be of the Fermi-Dirac form with two 
parameters (Ty,rjy), 



F2(7?,)r3exp(g/T,-7/,)-Kl' 



(12) 



where rji, is the degeneracy parameter, T^, is the neutrino 
temperature, and 



X dx 



exp(a; — 77) + 1 



(13) 



For numerical calculations, we will take {E^^) — 11 MeV, 
{E,^) = 16 MeV, {Ey^) = {E^^) - {E^^J = (E^^) = 25 
MeV, and r]^^ = ry^^ = rj^^ = 77^^ = 77^^ = jy^^ = 3. With 
these choices, we have T^^ ~ 2.76 MeV, Tp^ ~ 4.01 MeV, 
and = % = = Tp^ ~ 6.26 MeV. 

In principle, one could use the profiles of baryon den- 
sity nb, temperature T and electron fraction Ye (net 
number of electrons per baryon) obtained from numer- 
ical simulations of core collapse supernovae. Here we will 
use a simple analytical density profile, and approximate 
the envelope above the neutron star as a quasistatic con- 
figuration with a constant entropy per baryon S (see, 
6.17., Ref. [IBl)- Taking the enthalpy per baryon, TS, as 
roughly the gravitational binding energy of a baryon, one 
has the following temperature profile 



(14) 



'PI 



where Mns is the mass of the neutron star, ttijv is the 
mass of a nucleon, and mp\ ~ 1.221 x 10^^ MeV is the 
Planck mass. We assume that the entropy per baryon S 
in the hot bubble is dominated by relativistic degrees of 
freedom, 



2^ T3 



(15) 



where we have taken the Boltzmann constant fee and the 
reduced Planck constant h both to equal 1, and is 
the statistical weight in relativistic particles. Combining 
Eqs. (HH and ([T5]), one obtains the baryon density profile 



as 



27r2 
1^' 



S' 



'PI 



(4.2 X 10^" cm 



(16a) 



1.4M, 



100 



^lOkm^ ^ 



(16b) 



However, we note that, in reality, the baryon density 
rib near the neutrino sphere is much higher than that 
estimated from Eq. p6p . In fact, near the neutrino sphere 
the density profile is better represented by 



Tibo exp 



Ri, 



(17) 
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FIG. 2: (Color online) Plot of (effective) net lepton number 
density n; (r). The dashed and dot-dashed lines are for the net 
electron density rie — Yefib using the baryon density profile in 
Eq. Uni) with S = 140 and 250, respectively. The dotted line 
is for the net electron density assuming the baryon density 
profile in Eq. (|17p only. The solid line is for the efi'ective net 
Ve_ density along the radial trajectory [Eq. (|40|l ]. 



isothermal conditions in the vicinity of the neutron star 
surface. This, coupled with the expected very low elec- 
tron fraction Ye near the neutron star surface, impUes 
that the baryon density must have this exponential de- 
pendence on radius, at least for a radius interval ^ /ins- 
It turns out that addition of this exponential density 
profile near the neutrino sphere facilitates the multiangle 
simulations of neutrino flavor transformation. In Fig. [2] 
we plot the net electron number density 



(19) 



obtained from the exponential profile in Eq. p7|l . For 
comparison, we also plot r7,e(r) obtained from the con- 
stant entropy profile [Eq. p^ ] with entropy per baryon 
S = 140 and 250. In both Fig. [2] and in the rest of the 
paper, we take Mns = 1.4 Af©, = 11 km, Ye = 0.4, 
5s = 11/2, ribo = 1.63 x 10^^ cm^^ and /ins = 0.18 km. 
Note that once we have specified nbo and /ins our model 
for the physical environment in the hot bubble is com- 
pletely determined by the choice of entropy per baryon 
S. In units of Boltzmann constant per baryon, we expect 
S - 100 in the hot bubble [H. 



where nbo is the baryon density at the neutrino sphere, 
and 



tNS 



R 



' ?niatt(^i/) 



" V^/ns 



(0.052 km) 



R,, 



10 km 



1.4M 



Mns 



(18a) 



IMeV 



(18b) 



is the scale height with rinatt(-Riy) being the matter tem- 
perature at the neutrino sphere. This exponential fall-off 
in density is expected on general physical grounds and 
is found in, e.g., the Mayle and Wilson supernova sim- 
ulations 0. As discussed in Refs. [H, [2^, a steady 
state between neutrino heating and cooling results in near 



B. Neutrino Flavor Transformation in Supernovae 

Our objective is to study the flavor evolution of the 
neutrino field when and mix with neutrinos and 
antineutrinos of another active flavor (say Vr and Vr). 
We write the wave function of the flavor doublet of a 
neutrino (or antineutrino) as 



(20) 



where a and h are the amplitudes for a neutrino to be in 
the z/e {pe) and Vt (i?r) flavor states, respectively. The 
flavor evolution of ip^^ is determined by the Schrodingcr 
equation (see, e.g., Ref. Q) 



J 



-A COS 26* - 
A sin 29 



A + B A sin 261 + 



Acos29 - A - B 



(21) 



where 6 is the vacuum mixing angle. A, A and i3(er) are 
the potentials induced by neutrino mass difference, mat- 
ter, and background neutrinos, respectively. One obtains 
the appropriate Hamiltonian for antineutrinos by making 
the transformation 



A 



-A, 



B 



-B, 



B, 



-B* 



The vacuum potential is defined as 
A 



(23) 



where 6m is the neutrino mass-squared difference, and 
is the energy of the neutrino. (Note that we also 
use q as the energy or the magnitude of the momentum 
(22) of a neutrino in this section, which is the same as E^.) 
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We define the mass-squared difference in terms of tlie 
appropriate neutrino mass eigenvalues mi and 7713 to be 
Sm? = m| — nil. what follows we employ the nor- 
mal {Sm? = (Sm^^jjj) and inverted (Sm^ = — (5™^^^^) mass 
hierarchies. The matter potential is 

A = V2GFne = V2GFYenb, (24) 

where Gp is the Fermi coupling constant. We define a 
reduced density matrix g^, (in the flavor basis) from ip^, 



as 

^''-2\ 2a*b -|a|2 + |fe|2 J ' ^^^^ 

Note that this definition applies for neutrinos and an- 
tineutrinos. This is, however, different from the conven- 
tion adopted in Ref. Using Eq. (|25|) . the neutrino- 
neutrino forward scattering part of the Hamiltonian in 
Eq. (PT|) can be written as 



If B Ber 
2 I i?e% -B 



(26a) 
(26b) 



where q and q' are the momentum of the neutrino of interest and that of the background neutrino, respectively, and 
the flavor index is a = e or r. As mentioned above, neutrinos of the same initial flavor, energy and emission angle 
have identical flavor evolution. Consequently one must have 



(27) 



We note that 



J q • q'F(7?') dq' = J[smdsind\sin(j)sm(j)' + cos(j)cos(j)') + cosi}cos'd']F{i}')d{cosd')d(j)' (28a) 
= 27ry" cos?9cosi?'F(7?')d(cosi?'): (28b) 

where F(j!)) is an arbitrary function of i9, and we have used the cylindrical symmetry around the z-axis in deriving 
Eq. ((28b)) . Using Eqs. ©, (HID and (QH]), one can rewrite Eq. ((26b)) as 



H.. = ^^Y. /(I -cos 7? cos 7?') 



As noted in the introduction, previous simulations have used the single-angle approximation, wherein one assumes 
that the flavor evolution history of a neutrino is trajectory independent, 



((2II) 



and neutrinos on any trajectory transform in the same way as neutrinos propagating in the radial direction. Using 
the single-angle approximation, Eq. ((26b[|) can be further simplified to 



QyM')UM')j^^ - QljKq')UJyq')j^ 



dg', 



M M 



r 



where the geometric factor D{r/R^) is defined as 

2 



D{r/R,) EE i 



( — 

r 



(29) 



tion formalism (see, e.g., Ref. [l7l |) provides an intuitive 
way of understanding the results of our simulations. The 
wave function of a neutrino tpi, in Eq. ()20p can be mapped 
into a Neutrino Flavor Iso-Spin (NFIS) vector ^ using the 



Although our simulations are carried out by solving 
Eq. (I2ip numerically, the spin analogue of the wave func- 
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Pauli matrices cr: 



, , 2Re(a*6) 
- 2Ini(a*6) 
2 ^ |a|2-|5p 



1 



2Re(a6*) 
2Im(a6*) 

\a?-\b? 



(30a) 
(30b) 



((5m^ < 0). Because neutrinos can have different ener- 
gies, /xy varies from — oo to -|-cxd. The second term in 
Eq. ([33ll is induced by matter (neutrino-electron forward 
scattering), and we can write 



(36) 



and 



Note that the extra ay in Eq. (j30b[) transforms 2 of 
SU(2), the fundamental representation of antiparticles, 
into 2, the fundamental representation of particles. As a 
result, fo^ transforms in the same way as i;^^ under rota- 
tions. We also note the NFIS's defined in Eq. (|30p have 
constant magnitude 1/2. For a neutrino Va., ^v^z — +1/2 
(—1/2) for the pure z/g (I'r) state, where <;^^z is the 
third component of the NFIS. For an antineutrino Da, 
^vcz — +1/2 (—1/2) for the pure 9^ {v^) state. 

The NFIS <;^{q,'d) for cither a neutrino or an antineu- 
trino obeys the equation of motion 



2.^E/m^>^')w(.',^') 



(31) 



X /.'(?') 7ifTd(cosW 



where q and "d are the magnitude and polar angle of the 
momentum of the neutrino, H''^ is an effective field, and 
li{'d,'d') is the coupling coefficient between (;^{q,-d) and 
the background neutrino (;^i{q',i)') with 



fi{i},i}') = -2\/2G'f(1 - cos 1? COST?')- 



(32) 



The summation index v' in Eq. (j3ip runs over v,,, i't, 
and Vr_. According to Eq. (|3ip . the motion of a NFIS 
in fiavor space is analogous to that of a magnetic spin 
which simultaneously precesses around a "magnetic field" 



Me 



\/2Gf 



(37) 



Before we show the results of our simulations, we shall 
estimate the "MSW resonance radius" tmsw for a neu- 
trino with a typical energy. The MSW resonance condi- 
tion would be 



A cos 2Q — A(rMSw) 



(38) 



if we ignore the neutrino-neutrino flavor-diagonal poten- 
tial B. We will take \6m^\ = 3 x 10"^ eV^, the atmo- 
spheric value, and we will take the effective 2x2 vacuum 
mixing angle to be = 0.1. Note that this value is well 
below the experimental limit on ^13. For these parame- 
ters, the MSW resonance radius of a neutrino in the case 
of normal mass hierarchy [Srv? > 0) or an antineutrino 
in the case of inverted mass hierarchy ((5m^ < 0) with en- 
ergy = 10 MeV is tmsw - 127 and 59 km for S = 140 
and 250, respectively. We shall also estimate the radius 
for significant neutrino flavor transformation if neutrinos 
and antineutrinos are in the "synchronization" mode . 
When neutrinos are in the synchronization mode, all the 
NFIS's behave as one "magnetic spin" with 



2E, 



sync 



E 



, (39) 



where we have assumed all the NFIS's are aligned or an- 
H<=« and the other "spins". The "magnetic field" H'^** is tiahgned with e{. Because E^yr^c is positive, the sign 



composed of two components in our case, 
H'=ff(g) =/iv(g)Hv + MeHe 



(33) 



In Eq. p3p Hy stems from neutrino mass difference and 
can be written as 



Hv = -el sin 261 



e!, cos 29, 



(34) 



of the first term (left hand side) of Eq. (|39p should be 
chosen to be the same as that of the product 6m? {^y)- 
If (/iy) > and 5m? > 0, all the neutrinos and an- 
tineutrinos go through the same conversion process as 
a neutrino of energy E'sync at rMSw(-£'sync)- Similarly, 
if (/iy) > and 6m^ < 0, all the neutrinos and an- 
tineutrinos go through the same conversion process as an 



where e^^^^ are the orthogonal unit vectors in fiavor 
space corresponding to cr^(j,.2). Here we define 



antineutrino of energy E'sync at rMSw(£'s' 



Neutrino 



(35) 



where the plus sign is for neutrinos and the minus sign 
is for antineutrinos. With these definitions, neutrinos 
possess positive (negative) "magnetic moments" /xy and 
antineutrinos possess negative (positive) ones if Sm^ > 



fiavor transformation is suppressed for other synchroniza- 
tion scenarios. For the parameters we have chosen, we 
find that (^y) > if Sm'^ > 0, and (fiy) < if Sm'^ < 0. 
The characteristic energy of the synchronization mode 
is £'sync — 2.47 MeV for both cases. Therefore, in the 
synchronized mode neutrinos and antineutrinos should 



transform simultaneously at rMsw(^! 



sync J 



~ 80 and 37 

km for S ~ 140 and 250, respectively, if Sm^ > 0. These 
neutrinos/antineutrinos would experience very little fia- 
vor conversion if Sm'^ < 0. 
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A special case of synchronized behavior is the Back- 
ground Dominant Solution (BDS) JJl] where the NFIS's 
rotate in the plane spanned by and in flavor 
space. One of the necessary conditions for the BDS with 
large-scale simultaneous neutrino and antineutrino fla- 



vor transformation is that the flavor off-diagonal neutrino 
background potential B^^t dominates. To see this condi- 
tion more clearly, we define the effective net i^e number 
density along the radial trajectory as 



< = y (1 - 2 • q') dn,Jq')d<7' ^ J {1 ~ z ■ q') dn.M'W 
D{r/R,) 



(1.66 X 10-^'' cm"-^) 



1 - 



1-1^ 

r 



10 km 



LuJjlQ^^ erg/s) W(105i erg/s) " 
(i;^J/(10MeV) (£;pJ/(10MeV)_ 



(40a) 
(40b) 

(40c) 



Note that |Ber| \/2GF|n^f | and B = if the BDS 
obtains. (Obviously, if no flavor transformation has oc- 
curred, Ber = and B = V2Gvnf^.) We plot nj;^ (r) to- 
gether with ne{r) used in our simulations in Fig. [2l The 
neutrino background potential will dominate the mat- 
ter potential on the radial trajectory if > rie, which 
corresponds to a radius as low as r ~ 13 km for the 
parameters we have chosen. 

For numerical simplicity, we have fixed Yg = 0.4 in our 
simulations. Of course, the value of Ye actually varies 
with the radius and is affected by Ve and v,, fluxes through 
the weak interactions P, [2] 

+ n ^ P + e~ , (41a) 
De+p^n + e~^. (41b) 

The rates of these processes can also be affected by 
weak magnetism corrections [13, H^- In the nu- 
merical simulations presented below we have not in- 
cluded ncutrino/antineutrino flavor transformation feed- 
back through these processes on Yg- This is an important 
aspect of the physics of the supernova environment which 
we leave to a subsequent paper. 

Models for r-process nucleosynthesis can be sensitive 
to the value of in the region where T > 0.1 MeV 
For our chosen density profile, T = 0.1 MeV occurs at 
r ~ 139 and 78 km for entropy per baryon S' = 140 and 
250, respectively. These values of radius are well outside 
our simple estimates for where conventional MSW, syn- 
chronization, or BDS-like flavor conversion could occur. 
The numerical results to be discussed in the next sec- 
tion will give us a much better idea of where large-scale 
neutrino flavor transformation actually occurs. 



III. MULTIANGLE NUMERICAL 
SIMULATIONS 

In Sec. nil A1 we discuss our numerical calculations, and 
point out two potential pitfalls in any multiangle simu- 
lation. In Sec. nil Bl we show the main results from our 
multiangle simulations. For the simulation results pre- 
sented in this section we have taken |(5m^ | = 3 x 10~^ eV^, 
9 = 0.1, L^^ = Lp^ = L^^ = = 10^1 erg/s and 
S — 140 unless otherwise stated. 



A. Numerical Scheme 

We have developed two independent sets of numerical 
codes using different computer languages. We have used 
them to provide cross checks to obtain consistent results. 
Both codes employ a large multidimensional array of neu- 
trino wave functions ipi^c{Ev,cosdo) and V'Pq (-Ei/, cos??o) 
and evolve them simultaneously following the scheme 
outlined in Sec. |TT1 Each code employs an adaptive step 
size control mechanism, but the two codes have different 
ways of estimating errors and adjusting step sizes. The 
energy bins are chosen to have equal sizes for convenience 
in comparing neutrino energy spectra at different radii. 
The angle bins are determined in such a way that each 
bin has the same size in cost? at radius i?bin. In most 
cases we have taken i?bin — Rvi the neutrino sphere ra- 
dius. Note that the angle bins have different sizes in cos d 
\ir ^ i?bin- 

At the basic level, both codes obtain ■\\}u(l + 61) from 
'01^(0 by using the following equation as a first step: 
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+ SI) ~ exp{-iHSl)i^^{l) 

1 /Acos(A^/) — i/iii sin(A5/) 



A 



~ihl2 sm(XSl) 



— i/ii2 sin(A(5Z) 
Xcoa{XSl) + ift-ii sin(A(5/) 



Mi), 



(42a) 
(42b) 



r 




FIG. 3: (Color online) Average survival probability 
{PvcVe {^)) along the radial trajectory (cos i9o=l) with the nor- 
mal mass hierarchy in different numerical schemes. Here the 
average is done over the initial energy spectra of z/^. The dot- 
dashed line uses 160 angle bins and error tolerance 10~^ in 
each step without the initial baryon density profile n{^. The 
dashed and solid lines both include n(j and employ error tol- 
erance 10~^°, but use 256 and 512 angle bins, respectively. 
Calculations with 768, 1024 and 1407 angle bins in different 
binning schemes produce curves which fall on the solid line. 



where hn and hi2 are the diagonal and ofF-diagonal ele- 
ments of the Hamiltonian H . and 



A 



^ . /f,2 



hi. 



\h,2\'. 



(43) 



[Although not written out explicitly, the Hamiltonian H 
and its elements in the above equations have dependence 
on both the Afhne parameter I and trajectory angle i?, as 
can be inferred from Eqs. (pij) through (|26b[p .] We note 
that Eq. ([1^ preserves the unitarity of ipi, automatically. 
We also note that Eq. (|42)) becomes exact if H is inde- 
pendent of spatial coordinate. Therefore, the step sizes 
employed in our numerical codes are not restricted by the 
size of H but are restricted by the rate of change of H. 

In the course of our work we have discovered two pit- 
falls which apply to any multiangle scheme. Failure to 
avoid these pitfalls may lead to quantitatively or quali- 
tatively inaccurate results (see Fig. [3]). 

The first potential problem has to do with the expo- 
nential term in the profile for the baryon density [see 
Eq. (|17p ]. The baryon density is very high near the neu- 
trino sphere when n'-^ is included. This sometimes forces 



numerical schemes to employ initially very small step 
sizes. The numerical codes using the single-angle approx- 
imation can generally drop n'-^ without loss of accuracy 
at large radius. These codes will of course run faster 
without n'^. However, in multiangle simulations, ignor- 
ing n(j makes the background neutrino potential B much 
bigger than the matter potential A even at the neutrino 
sphere. As a result, the evolution histories of neutrino fla- 
vors on all trajectories are strongly coupled starting from 
the beginning. This strong correlation among all trajec- 
tories also forces small step sizes. In addition, without 
n(j there is a tendency for neutrinos to undergo flavor 
transformation very close to the neutrino sphere. This 
behavior is suppressed if there is a large and dominant 
matter potential A. Including makes the matter po- 
tential A(Ri,) much bigger and this helps keep neutrinos 
in their initial flavor states, at least for the significant 
range of neutrino/ antincutrino energies and for our cho- 
sen value of l^m^l ~ Sm'^^^^^. We also note that neutrinos 
on different trajectories propagate through different dis- 
tances. A big matter potential breaks the correlation 
between neutrinos on different trajectories and lets them 
evolve independently for awhile. 

These considerations can be cast in simpler, more phys- 
ical terms. In the relatively narrow region near the 
neutrino sphere where n'^ dominates it has the effect of 
changing, or "resetting" , the neutrino wavefunctions rel- 
ative to what they would have been had we employed the 
unphysical low-density profile all the way to the neutrino 
sphere. In the latter unphysical case, neutrinos are in fia- 
vor eigen states at the neutrino sphere, and the NFIS's 
are perfectly aligned with each other yet slightly devi- 
ated from the total effective field H'^*^. The effects of 
this unphysical setup does not go away quickly with in- 
creasing radius because the coupling among the NFIS's 
(arising from neutrino-neutrino forward scatterings) is so 
strong. If the exponential baryon density profile n'-^ is 
added, the overwhelming matter field He at the neutrino 
sphere not only makes the NFIS's more aligned with H°*^, 
but also breaks the coupling of the NFIS's propagating 
along different trajectories. In the short distance where 
the matter field He dominates, the NFIS's on different 
trajectories have traveled different distances and so have 
developed different phases. At the radius where n'^ be- 
comes negligible, the NFIS's are effectively "reset" to a 
more physical condition than one would obtain without 

<■ 

The other pitfall is that one may use an insufficient 
number of angle bins. Assuming that there has been 
very little neutrino flavor conversion close to the neutrino 
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sphere where r ~ R^,, we can write 



B{r,d) 



2^Rl \{E,^ 
Be.(r,7?) ~0. 
For a small step size SI, one has 

^^Jl + SI) ~ cxp(-iiJ(5/)V'^„ (/) 




-i^^sin[(A + B)5/]^ 

i{A+B)Sl 



(44a) 
(44b) 

(45a) 
(45b) 



where we have used the fact that A + B ^ A at r ^ R,^. 
It is the off-diagonal elements of the transformation ma- 
trix in Eq. (|45bl) that govern the exchange of the two 
flavor components of a neutrino wave function. These off- 
diagonal terms, we note, are oscillatory functions of the 
B potential and the step size SI, both of which have angu- 
lar dependence [see Eqs. (0) and (|44ap ]. Physically, this 
oscillatory feature with respect to angles is suppressed by 
strong correlation among neutrinos on different trajecto- 
ries. Numerical codes without enough angular resolution, 
however, could allow a spurious "cross-talk" between an- 
gle zones which artificially strengthens flavor oscillations. 
This unphysical feedback could produce substantial neu- 
trino flavor conversion even at low radius in some numer- 
ical schemes. 

In Fig.[3]we plot average survival probability {Py^u^ {f)) 
along the radial trajectory with the normal mass hierar- 
chy using different numerical schemes (error tolerance, 
number of angle bins, etc.). Here P^^^^(r) is the prob- 
ability for a Ve to be a i/g at radius r, and the average 
is done over the initial energy distribution for Vf,. (As 
mentioned above, we use Va and to denote the neu- 
trinos and antineutrinos that are emitted in flavor state 
a at the neutrino sphere.) One sees that spurious neu- 
trino flavor transformation (dot-dashed line) could occur 
at low radius with a combination of insufficient number 
of angle bins, a loose error control, and neglect of n'^. 
If we employ = 10^-'^ erg/s and choose a stringent 
error tolerance (~ 10"^") at each step, we find that it 
takes > 500 angle bins in order to achieve convergence 
and run-to-run consistency. Because the B potential in- 
creases with neutrino luminosity, we expect that even 
more angle bins would be required to obtain convergence 
at larger neutrino luminosity. 

Our numerical simulations generally employ > 500 an- 
gle bins and > 500 energy bins for each neutrino species. 
Typically, our codes execute > 10^ steps during each pro- 
duction run. It is clear that multiangle simulations are 
only feasible using large-scale parallel computation. 



B. Simulation Results 

In Fig. HJa), we plot (-Pi/^i/^ (r)) with the normal neu- 
trino mass hierarchy (Sm? > 0) on both the radial 
(cos7?o = 1) and tangential (cost?o = 0) trajectories. For 
comparison, we also plot {Pu^uA''')) for the = Q {A po- 
tential only) case, which is obtained from the single-angle 
simulation by setting = 0. The = case corre- 
sponds to the limit where neutrinos go through MSW 
resonances independently of each other. In the full syn- 
chronization limit, all neutrinos and antineutrinos un- 
dergo flavor transformation in the same way as does a Uf. 
with energy -Esync in the standard MSW mechanism. Us- 
ing only the matter potential, we have calculated P^^^,^ 
for a Ve with energy -Egync propagating along the radial 
trajectory. The result is shown in Fig. HJa). The results 
of our simulations are clearly different from those in the 
= and full synchronization limits. In particular, 
our simulation has {Pv^u^ir)) crossing the 1/2 line later 
than in the synchronization case, but earlier than in the 
Li, = case. We also note that {Pi^^iy^{r)) oscillates and 
even bounces back after an initial decrease. 

In Fig.[4ljb) we plot (Pu^p^ir)) in the normal neutrino 
mass hierarchy scenario. The conventional MSW conver- 
sion of De is suppressed if Sm"^ > 0, which is illustrated 
dramatically by the L,y — case. In the case of full 
synchronization, will be converted into Dr in exactly 
the same way as Ve is converted into i^r- The results 
of our simulations are again like neither of these limits. 
Unlike the L,y — case, the actual values of {Pv^v^ir)) 
may substantially decrease at some values of radius, and 
unlike the full synchronization case, {Pu^uA''')) oscillates 
and bounces back to nearly unity at large radius. 

The results of the inverted neutrino mass hierarchy 
[Sm? < 0) are more surprising. These are plotted in 
panels (c) and (d) of Fig. [4l The full synchronization 
limit predicts no flavor conversion for both and Pg, and 
the = limit predicts that only antineutrinos will be 
converted. Our simulation finds substantial conversion 
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FIG. 4: (Color online) Plots of average survival probability {P^^v^) (left panels) and {Pv^a^) (right panels) with the normal 
(upper panels) and inverted (lower panels) neutrino mass hierarchies, respectively. The solid and dotted lines give average 
survival probabilities along trajectories with cosi^o = 1 and cosi9o = 0, respectively, as computed in the multiangle simulations. 
The dot-dashed lines and the dashed lines characterize the limits where neutrinos and antineutrino undergo flavor transformation 
individually (L^ = and A potential only) and simultaneously (full synchronization), respectively. The dashed line is not 
distinguishable from the dot-dashed line in panel (c). 



of both Vf, and v^. Furthermore, this phenomenon occurs 
at a radius even smaller than that expected in the full 
synchronization limit with 5m? > 0. Again, we note that 
in the inverted mass hierarchy scenario both {Pi,^v^{r)) 
and {Pp^iy^{r)) oscillate after flavor transformation starts. 

In Fig. [5] we plot P,y^i,^{r) for ly^ with a few character- 
istic energies on both the radial and tangential trajec- 
tories. We have employed the normal mass hierarchy in 
this calculation. One sees that the Piy^i,^{r) curves have 
similar trends with radius over most of the Ve energy 
range considered. This is especially true for the values 
of radius where neutrino flavor transformation has just 
become significant and for the tangential trajectory. 

The results presented in Fig. [5] lead us to conclude 
that the flavor transformation histories of neutrinos on 
different trajectories can be very different. To illustrate 
this point more clearly, we plot in Fig.ini^a) P^^^^{cos-do) 



at r ~ 92.85 km for j/g with specified energies, and 
employing the normal mass hierarchy. Indeed the val- 
ues of Py^^^ (cos i^o) vary with angle, especially around 
cos 1^0 — 1- Moreover, the trend of Pi,^^^ {cos -do) with 
angle is similar for i^e with different energies over most 
of the energy range considered. This again demonstrates 
the collective feature of the neutrino flavor transforma- 
tion in the hot bubble. 

In Fig. El^b) we plot the corresponding antineutrino 
survival probability Pp^p^ (cos zJq)- This also shows angu- 
lar dependence and collective flavor transformation. In 
Fig. Eljc-d), we plot Pi/^jy^ (cos t9o) and Pp^p^ (cos i9o) with 
the same parameters as in panels (a-b) but at r ~ 88.57 
km and with the inverted mass hierarchy. It is interest- 
ing to see that, in addition to the features pointed out 
for panels (a-b), in the inverted mass hierarchy case both 
Pp_. (cos "iJo ) and Pp_.p^ (cos i^o) oscillate over most of the 
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FIG. 5: (Color online) Plots of P^^,^^(r) with the normal mass hierarchy. Panel (a) is for the radial trajectory (cost?o = 1), 
and (b) is for the tangential trajectory (cosj9o ~ 0). The dot-dashed, dotted, dashed and solid lines are for Ve of energies 6.95, 
8.95, 10.95 and 14.95 MeV, respectively. 




FIG. 6; (Color online) Plots of P^^i^^ (cos i9o) (left panels) and Pv^p^icos'&o) (right panels) at r ~ 92.85 km for the normal mass 
hierarchy (upper panels) and at r ~ 88.57 km for the inverted mass hierarchy (lower panels). The dot-dashed, dotted, dashed 
and solid lines are for Ve or Ue of energies 6.95, 8.95, 10.95 and 14.95 MeV, respectively. 
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FIG. 7: (Color online) Change of energy spectra of neutrinos (left panels) and antineutrinos (right panels) with the normal 
(upper panels) and inverted (lower panels) neutrino mass hierarchies. The dotted and dot-dashed lines are the spectra of 
neutrinos (antineutrinos) in the electron and tan flavors, respectively, at r = i?^, and the solid and dashed lines are the 
corresponding spectra at r = 250 km. 



range of cos??o- 

In these simulations, significant neutrino flavor trans- 
formation ends at r '--^ 230 km (Fig. [4]). To see how the 
energy spectra of neutrinos and antineutrinos have been 
altered by flavor transformation, in Fig. [7Ja) we plot both 
fvS^) ^^'^ ft^-ri^) the neutrino sphere, and fv^{E) 
and (E) a,t r — 250 km. Here we have employed the 
normal mass hierarchy and we take fu{E) to be propor- 
tional to both fi^{E) and the flux of v [e.g., fvS^) ^ 

UAE) Ea / P'^^^An.^ and U^E) « UAE)lJ{E,J], 
such that 

Here the scheme for angle- aver aging the energy spectra is 
simply the angle dependence in the neutrino flux "seen" 



by a nucleon at radius r. As a result, the angle-averaged 
spectra shown are those appropriate for use in the weak 
interaction rates. It is interesting to see that most of 
the low-energy {Ei, < 9.5 MeV) i^e are converted into 
i^r, while a significant fraction of high energy survive. 
We also plot the corresponding energy spectra of J>e and 
D-r in Fig. [TJ^b) . The energy spectra of antineutrinos are 
changed very little in the normal mass hierarchy scenario. 
The energy spectra of neutrinos and antineutrinos in the 
inverted mass hierarchy scenario are plotted in Fig. [7fc) 
and (d), respectively. In these figures, both and De 
swap spectra with Vr and i>T- , respectively, over a signifi- 
cant energy range. 



The numerical results that we have presented cannot 
be explained easily by the conventional MSW mechanism 
or by synchronization. We will try to develop some in- 
sight into, and understanding of these results in the fol- 
lowing section. 
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IV. SINGLE-ANGLE SIMULATIONS AND 
PHENOMENOLOGICAL ANALYSIS 

To understand the numerical results obtained from the 
multiangle simulations, we have re-examined the numeri- 
cal simulations using the single-angle approximation with 
similar setups and initial conditions. We found that al- 
most all the interesting features seen in the multiangle 
simulations are also present in the single-angle simula- 
tions, though they can differ in a quantitative sense. The 
simulations performed using the single-angle approxima- 
tion do not have the numerical difficulties that are the 
hallmark of the multiangle ones, and they require fewer 
computational resources. Most importantly, the single- 
angle simulations produce results qualitatively similar to 
those in the multiangle simulations, and yet do not in- 
volve complicated entanglement of neutrino flavor trans- 
formation on different trajectories. They are therefore 
easier to understand. In Sec. IIV Al we will try to ex- 
plain some of the results presented in Sec. IIIIBI with the 
help of these simplified calculations. In Sec. IIV B] we will 
study how the onset of large-scale collective neutrino fla- 
vor transformation is related to the neutrino luminosity 
Li, . We will comment on the validity of the single-angle 
approximation at the end of this section. 

Unless otherwise stated, all the simulations discussed 
in this section have the same parameters as those in 
Sec. mi i.e., \5rr?\ = 3 x 10-3eV^ 9 = 0.1, = lO^i 
erg/s and S = 140, but are based on the single-angle 
approximat ion . 



direction. Because the densities of neutrinos and antineu- 
trinos are also large in this region, the i?er potential in 
Eq. (|2ip dominates, and both neutrinos and antineutri- 
nos are in a state similar to the Background Dominant 
Solution (BDS) TG]. 

The numerical results clearly have shown that neu- 
trinos and antincutrinos undergo some collective flavor 
transformation in the hot bubble with the neutrino mix- 
ing parameters we have used. The collective modes of 
flavor transformation that neutrinos may have in the hot 
bubble, according to Ref. [13], are either the synchro- 
nization or the bi-polar type. FIGs. [4] and [8] show that 
the collective mode corresponding to the conditions and 
parameters used here does not conform to the full syn- 
chronization limit. Therefore, we will focus our discus- 
sion on the bi-polar flavor transformation. Neutrinos and 
antincutrinos can have substantial flavor transformation 
simultaneously only through the bi-polar mode in the in- 
verted mass hierarchy scenario. This is seen in our simu- 
lations. In general, the region where neutrinos transform 
through the bi-polar mode is characterized by parameters 
satisfying [l3| 



< K < 



(48) 



where e is a measure of the difference in the energy dis- 
tribution functions of i>e + i^r and De + i^r, the parameter 



^ \Sm^\/2{E,) 



(49) 



A. Neutrino Flavor Transformation in the Bi-Polar 

Mode 

The novel features of neutrino flavor transformation in 
the hot bubble region are easier to understand in the for- 
malism of NFIS (Neutrino Flavor Iso-Spin) [l3] than in 
the traditional formalism of the wave functions. In Fig. [SI 
we plot (<^3;(r)), {<;y{r)) and (^^(r)), the three components 
of the average NFIS's in flavor space, for and in both 
the scenarios with a normal mass hierarchy and with an 
inverted mass hierarchy. (The three components of the 
NFIS's are averaged over the initial neutrino or antineu- 
trino energy spectra.) We note that the probability for a 
neutrino or antineutrino initially in the a flavor state to 
be in the electron flavor state is related to <,z by 



P- _ 



1 



2 



(47a) 
(47b) 



Comparing Fig. [8] with Fig. [H one sees that the results 
of single-angle simulations are qualitatively the same as 
those obtained in the full multiangle simulations. We 
also note that in the region where neutrinos transform, 
the NFIS's of both neutrinos and antincutrinos have large 
values of and c^y, and roughly precess around the 



gives the strength of background neutrino effect through 
the effective single species neutrino number density 
[e.g., nf{L^,r) ~ D{r / R^){L^/ {E^)) /2ttRI for the ra- 
dial trajectory], and 5E^ is the characteristic width of 
the neutrino energy distribution. 

Neutrinos are in the synchronization mode if e > k. 
This corresponds to the high neutrino luminosity limit. 
Neutrinos will transform individually through the MSW 
mechanism if k > {E^) /25E^, which is effectively the low 
neutrino luminosity limit. Using the single- angle approx- 
imation, one has [171 ] 
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For simplicity, we assume that {E^^) and {E^^) charac- 
terize the energies of v,, -\- and + i^n respectively, 
and [13 



{{E^J^{E,jr 
2{{E^J^ + {E,J^) 



0.033. 



(51) 



15 




T 1 1 1 r- 




II ill 



(a) 



_l I I I I 1_ 



100 



_i I I I i_ 



_i I I I I I i_ 



0.4 



0.2 



-0.2 



-0.4 



— 1 1 1 1 1 1 1 r- 



200 



Ml I 

Ml ^ , 

■ I MM:[i^^^^il!--WV^•■■ 
l:|i|; i^-il^-^- 



200 




(b) 



100 



r (km) 



r (km) 




FIG. 8: (Color online) Plots of {<;vcir)) (left panels) and {'ii?s(r)) (right panels) with the normal (upper panels) and inverted 
(lower panels) mass hierarchies, respectively. The dashed, dotted and solid lines represent {<;x{r)), {<;y{r)) and (?z(r)), respec- 
tively. Note that {Pv^u^) = 1/2 + (c^^z) and {Pa^v^) = 1/2 — (Ci^^z). The insets in panels (c) and (d) are blowups of the 
corresponding plots in the range 62 km < r < 66 km. These are single-angle calculation results. 



Using Eqs. ([35]), (|50bp and ([5T|) . we estimate that neu- 
trino flavor transformation exits the synchronization 
mode and enters the bi-polar mode (Bi-polar Starting) 
at TBS ^ 73 km for the parameters we have used. Taking 
5E^/{Ey) ~ Tyj^Ey) ~ l/4, we estimate that collective 
neutrino flavor transformation ends (Bi-polar Ending) at 
tbe ~ 202 km. Beyond this point conventional MSW 
flavor transformation takes over. 

Although these values are crude estimates based on 
oversimplified assumptions, we find that they roughly 
match the region where neutrinos and antineutrinos 
transform simuhaneously in the inverted mass hierarchy 
scenario. Therefore we conclude that collective neutrino 
flavor transformation observed in our fuU numerical sim- 
ulations is indeed of the bi-polar type, as predicted in 
Ref. [H]. 

An interesting feature of neutrino flavor transforma- 
tion in the bi-polar mode is that the transformation is 



not completely suppressed by the large matter potential 
A in Eq. ([21]) if < 0. From Fig.|4]and Fig.[8l one sees 
that the flavor transformation with 6m^ < may occur 
at values of radius even smaller than those predicted in 
the synchronization limit with 6m'^ > 0. Although collec- 
tive neutrino flavor transformation in the bi-polar mode 
has been studied in the zero and large matter potential 
limits [iTj . an analytical or semianalytical analysis has 
yet to be performed to show how neutrinos transform in 
the bi-polar mode as the matter potential A decreases 
and approaches the vacuum potential A. 

In Fig. [UJa) we plot i^u^ziE^^) in the normal mass hi- 
erarchy scenario at 400 km for the cases with L^, = 10^^ 
and 5 x 10^^ erg/s. One immediately sees that there is a 
rather sharp transition edge at Eq ~ 9.5 and 7.9 MeV for 
Ly = 10^^ and 5 x 10^^ erg/s, respectively. Noting that 
Pv^u, = l/2 + c:,y^z, one sees that with energy below Eq 
are almost fully converted into v^, while with energy 
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FIG. 9; (Color online) Plots of <;,^^z{Ei,^) (left panels) and <;i?^z(Ei^^) (right panels) for both the normal (upper panels) and 
inverted (lower panels) mass hierarchies at r = 400 km. The dashed and solid lines are for = 10^^ and 5 x 10^^ erg/s, 
respectively. These are single-angle calculation results. 



above Eq but below another threshold En mostly sur- 
vive. The threshold Eh is roughly at ~ 22 and 40 MeV 
for Li, — 10^^ and 5 x 10^"'^ erg/s, respectively. Because 
Pv^u^{E) — Py^i,^{E), there exists a similar transition 
edge for Vt- This difference in the flavor transformation 
of the neutrinos of low and high energies is responsible 
for the partial swap of the spectra of Uf. and v-r seen 
in Fig. Elja). We also plot the corresponding values of 
^i?^z{Ev^ in Fig.[9^b). Knowing that Pp^p^ = 1/2 - fe^^, 
one sees that most of the survive. This is also true for 

We plot <ii,^z{Eu^) and (^p^ziEp^) with the inverted mass 
hierarchy in Fig.[9jc) and (d), respectively. There also we 
see a transition edge with Eq — 8.5 MeV for v^, which 
is similar to that with the normal mass hierarchy, but 
reversed in direction. We note that Eq is essentially the 
same for both = 10^^ and 5 x 10^^ erg/s in the case of 
an inverted mass hierarchy. This transition edge results 
in the partial swap of the spectra of Ve and Vr shown 
in Fig. ^c). The behavior of De is more complicated. 



Roughly speaking, with energy below some threshold 
i?L or between Em and En are mostly converted into Vt , 
where =i 3 and 1.8 MeV, Em ^ 16.5 and 8 MeV, and 
Eh ~ 20 and 40 MeV for U = 10^^^ and 5 x 10^^ erg/s, 
respectively. 

As mentioned above, we do not have an analytical or 
semianalytical analysis of the bi-polar mode flavor trans- 
formation in the general cases. Nevertheless, we propose 
a tentative explanation of the main features in Fig. [5] as 
follows. 

The results shown in Fig. [8] suggest that the NFIS's of 
both neutrinos and antineutrinos roughly rotate around 
§2 with a frequency uj. In fact, in the limit that A A, 
the NFIS's in the bi-polar mode rotate around the vac- 
uum field Hv = = sin 261 -He^ cos 261 [see Eq. ([M)) ]. 
which is close to if ^ 1 . 

To give a rough feel for the behavior of such a system, 
let us study a toy scenario where ^(i?,/) is coupled to both 
Hy and another field which rotates in the plane 

perpendicular to Hy. Thus the equation of motion for 
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FIG. 10; (Color online) The toy scenario explaining the evolution of NFIS <f. The NFIS <; can be viewed as a "magnetic spin" 
which is coupled with a constant field Hv and a field S rotating with angular frequency ui. We actually consider here two 
"magnetic spin" couplings, one with "magnetic moment" fiv and the other with jj,^. The vector /x^S rotates in the e^-e^ plane 
in the clockwise sense when viewed from above, looking in the — direction. The "magnetic spin" <j is aligned initially with 
the dominant field S at time t = [panels (a) and (d)]. The problem is easily solved in the corotating frame where S is not 

rotating. In this corotating frame, the "magnetic spin" <j rotates around the effective field H = {fiv — u>)ilv + [panels (b) 

~ efF 

and (e)]. If S slowly reduces its length to zero, the angle between the spin <f and the effective field H is constant (adiabatic 
process) , and <; ends up almost aligned with Hv at i = co if — c^; > [panel (c)] or almost antialigned with Hv if /iv —uj<0 
[panel (f)]. 



<i{E^) can be written as 

+ /i,yS(i)(eJ^ coswt — smcji)], 

where /i^ is a coefficient, and e^, and are a set 
of orthogonal unit vectors with e^^ ^ % = ^z- Suppose 
that at i 0, we have \n^Y.{t = 0)| > |^v(£^i.)Hv| 
and <;{E,y) is aligned or antialigned with S(i — 0). We 
want to find the configuration of <i{E^) as E(i) slowly 
decreases toward zero. Eq. ([5^ turns out to be very 
simple in a corotating frame in which S(t) is fixed in one 
direction, say w. (Ref. [l3| points out the utihty of the 
corotating frame.) The equation of motion of <i{E^) in 
this corotating frame is 

^AE) = ^{E,) X H^ff (53a) 
at 

= ^{E,) X [ME,) - c^)Hv + /^.S(i)w], 

(53b) 

where a vector with a tilde symbol is the same as that 
without but viewed in the corotating frame. As E(t) 
decreases, rotates from the direction of /i,yW to that 
of (/xv (-£-!/) ^^)Hv. If this process is slow enough, ^{E,y) 
stays aligned or antiahgned with H°^, depending on the 
initial conditions, and will be either ahgned or antialigned 



with Hv when E approaches zero. We define 

w{E,) = ^i,[^iy{E,) - uME,) ■ S]t=o- (54) 

One can check that ^{E,), and therefore <;{Ei^), will be 
aligned with Hy as i —> cxd if > 0, and will be an- 
tialigned with Hy if w < 0. There can be a sharp tran- 
sition in the orientation of <; at energy E, = Eq, where 
/iv(i?c) — ^- The general features of this toy problem 
are shown in Fig. [TUl 

This analysis applies to collective neutrino flavor trans- 
formation in the hot bubble if (1) neutrinos and an- 
tineutrinos are in the collective mode even in the region 
where ^ <C A, and (2) the frequency of rotating NFIS's 
varies significantly more slowly than the neutrino density 
n,. In this case, S(i) corresponds to the rotating total 
NFIS, which decays as the neutrino density goes down 
with increasing radius. Because dominates in num- 
ber over other neutrinos and antineutrinos, the factor 
[<i{E„) ■ S](=o in Eq. ([51)1 is essentially the scalar product 
of the NFIS of the neutrino in question and that of the 
total ^1/^, which is positive for Ve and negative for De- For 
the normal mass hierarchy {5m? > 0), one has w > 
[note this behavior in Fig. EJa) and (b)]. Noticing that 
/i^ < [Eq. ([321)], one finds that ru is always negative for 
fee [Eu^ ) , which will be antialigned with Hy ~ in the 
end, as we have seen in Fig. [5] (b). One has tu < for 
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<iv^ {E^J if E,y^ < Eq and tn > if E^^ > Ec, where 



Er 



5m? 



2uj 



(55) 



condition 



'«(»'be) 



We see that is either approximately —1/2 or +1/2, 
depending on whether E^^ is less than or greater than 
Eq- This behavior can be seen in Fig. [9] (a). For the 
inverted mass hierarchy {5m? < 0), one has w < [see 
the small insets in Fig.[8l^c) and (d)]. One finds that w is 
always positive for (E'p^ ) , and will be roughly aligned 
with in the end, as we see in Fig. [9] (d). For <;^^ {E^^), 
one has w > if Ejj^ < Eq, and tu < if i?^^ > Eq- 
Therefore, the corresponding (^u^ziE^^) transitions from 
+1/2 to —1/2 as Ei,^ increases and crosses Ec, as we see 
in Fig. El (c). 

The above reasoning is, however, based on an idealized 
case. In reality, some NFIS's of high energy may never be 
locked into a collective bi-polar mode with other NFIS's 
under some conditions. Some NFIS's of moderate en- 
ergy may start to peel away from the bi-polar mode in 
the region where the matter potential A is comparable 
to A. In addition, some NFIS's of high energy may go 
through the conventional MSW conversion after the col- 
lective mode breaks down. Our guess is that Ve and De 
of energy E^, > En in Fig. [5] (a) and (d) have A > A 
when the collective mode breaks down, and they are at 
least partially converted through the MSW mechanism. 
The i>e with E^, < Ei^ in Fig. [5] (d) never enter the bi- 
polar mode, and are converted to Vt through the MSW 
or synchronization mechanisms. The i>f, with energies be- 
tween E'l and Eyi may have complicated flavor evolution 
histories which quite early cease to follow the bi-polar 
mode. 

Our argument becomes more accurate at high neutrino 
luminosity. With larger L^, more low-energy neutrinos 
and antineutrinos join the bi-polar flavor transformation, 
and more of them are locked into this collective mode 
until A < A. As a result, the threshold energies i?L and 
Em decrease, and £'h increases as Ly goes up. This is 
indeed the case as one can see from the comparison of the 
simulations with L^, — 10^^ and 5 x 10^^ erg/s (Fig. (2]). 

We have assumed to be a constant in our idealized 
analysis. This is not the case in reality. From Fig. [5] 
one sees that \lo\ slowly decreases with radius. If Ly is 
large enough, neutrinos and antineutrinos will be in the 
bi-polar mode even at values of radius where the mat- 
ter potential A{r) is negligible. We expect to be a 
function of (5m^, 9, f^{E^) and the local neutrino density 
n'^^{Li,, r), but to be independent of S, Y^, etc.. We note 
that neutrinos and antineutrinos start to deviate from the 
collective mode behavior at some radius rc as^ <f adiabat- 
ically rotates away from the direction of /XjyS. Further, 
we note that the value of Eq should be determined from 
uj(rc) using Eq. (j55p . One can attempt to estimate 
(< tbe) from Eqs. (pS)) and (|49p directly, resulting in the 



l^m^l 1 (E,) 

4V2GFnf{L^,rBE){Eu) ~ 26E^ 



(56) 



The value of tbe derived from Eq. (|56|) is an overesti- 
mate of rc- We have seen in Fig. [5] that collective flavor 
transformation ceases at r > tbe- However, at r ~ rc, 
all the NFIS's begin to slightly deviate from alignment, 
but are more or less still following the collective mode. 
Nevertheless, we expect K(rc), like K(rBE), to be deter- 
mined by fv{E^) only. As a result, nl^{L^, rc), and thus 
(jj(rc) and Ec, are actually independent of L^, if is 
large enough. 

We have calculated the energy spectra of neutrinos at 
r 400 km using the single-angle approximation with 
S = 140 and 250, and = 10^^ and 5 x lO^^ erg/s. The 
values of Ec in most of the cases agree well with each 
other for the same neutrino mass hierarchy. The value 
of Ec in the case with S = 140 and — 10^^ erg/s 
is different from those in the other three cases for the 
normal mass hierarchy [see, e.g., Fig. [Hla)] because L^, 
is not large enough, or equivalently, the baryon density 
profile is not sufficiently condensed toward the surface of 
the neutron star. We also note that ^^^z{Eu^) is not a 
strict step function, but has a transition region of finite 
width. The transition region in the normal mass hier- 
archy scenario overlaps with that in the inverted mass 
hierarchy scenario, which seems to suggest that the val- 
ues of |a;(rc)| are at least similar in these two cases. 



B. Onset of Collective Neutrino Flavor 
Transformation 

The radius where significant neutrino flavor transfor- 
mation starts can be very important for nucleosynthesis 
and for estimates of the expected late-time neutrino sig- 
nal [2ii[33, We define rx as the radius where {Pu^^J) 
falls just below 0.9. In Fig. [TITa). we plot rx{L^) for the 
cases with S — 140 and 250 in the normal mass hierarchy 
scenario based on our single-angle and multiangle simu- 
lations. For both entropy values, rx{L^) in single-angle 
simulations monotonically decreases as Li, increases. As 
a comparison, we also plot the corresponding values of 
rx/MSw{Esync) in Fig. [IHa). Here E,y^, ~ 2.47 MeV 
is the characteristic neutrino energy for the full synchro- 
nization mode, and rx/MSw(^i^) is the radius where a 
with energy Ei, has P^^u^ = 0.9 in the standard MSW 
mechanism. One sees that the values of rx{Lu) asymp- 
totically approach rx/]y[gvv(^sync)- This is not a surprise. 
According to Eq. ((48|) . neutrinos are in the synchroniza- 
tion mode if n'^J^{Ly,r) is large. In turn, n^^{Li^,r) in- 
creases with increasing L^, at a fixed radius r. As L^, 
increases, more and more low-energy neutrinos and an- 
tineutrinos are locked into the synchronization mode, and 
the characteristic neutrino energy of the synchronization 
mode decreases and asymptotically approaches -Egync- 
One also sees that for the same L^, the radius rx{L^) is 
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FIG. 11: (Color online) Plots of ry^{Lu) for both the normal [panel (a)] and inverted [panel (b)] neutrino mass hierarchies. 
The dashed and solid lines are based on single-angle simulations with S = 140 and 250, respectively. The cross and star 
symbols are based on the average of rx on the radial and tangential trajectories in our multiangle simulations with S — 140 
and 250, respectively. The error bars associated with the cross and star symbols indicate the range of values of rx on 
different trajectories. These are too small to be visible in most of the cases. For the normal mass hierarchy case [panel (a)], 
rx asymptotically approaches rx/MSw(^'sync) with large L^, which is 73.9 (dot-dashed lines) and 34.4 km (dotted lines) for 
S — 140 and 250, respectively. The dot-dashed line in panel (b) represents a crude estimate for where the collective neutrino 
flavor transformation exits the synchronization mode and enters the bi-polar mode [Eq. (|57[) ]. 



much closer to ?'x/MSw(^sync) in the S = 250 case than 
in the case virith S = 140. This is because with larger 
S the baryon density profile is more condensed toward 
the neutrino sphere, and ry^{Ly), hke rx/MSw(-E'sync), is 
smaller in this case. Therefore, n°^{L^, rx) is larger with 
a larger S but the same L^, and the synchronization is 
more complete. 

It is interesting to see that the values of rx{Lu) 
obtained from the multiangle simulations all fall be- 
tween those from the single-angle simulations and 
''X/Msw(£'sync)- Comparing Eq. (|26b[|) with \2Qhn . 
one can see that the single-angle approximation uses 
rfJ^{L^,r) on the radial trajectory. Note that nl^{L^,r) 
has smaller values on the radial trajectory than it does 
on any other trajectory. On average, values of rfJ^{Ly, r) 
are larger in the multiangle simulations than in the single- 
angle ones at the same radius r. At the same time, the 
full synchronization mode obtains when n1^{L^,r) 
oo. Thus rx{L^) computed from single-angle calcula- 
tions gives upper bounds on the actual rx{L^), and 
''x/MSw(-^sync) givcs a lowcr bound on this quantity. 

In Fig. IllT b). we plot the numerical values of rx{Lv) 
in our single-angle and multiangle simulations for the 
cases with S = 140 and 250, respectively, and employing 
the inverted mass hierarchy. One sees that the values 
of rx(ij/) monotonically increase with L^. In addition, 
they are not very sensitive to the value of S. To explain 
this phenomenon, we note that res , the radius where the 
neutrinos exit the synchronization mode and enter the 



bi-polar mode, can be estimated from the condition 

4V2Gf nf{L,,T^s){E.) " 2{{E,^)^ + {E,J^) 

[see Eqs. (08]), (|49]) and ^]. Clearly nf{L„, tbs) de- 
pends only on Sm^ , 6 and f,y{E^). Once these param- 
eters are specified, n'^^ {L^,rBs) is fixed. As a result, 
fBsiLi/) must increase with L^, for a fixed neutrino den- 
sity nl^{L,y,rBs)- In addition, rBs{Lv) depends only on 
Sm?, 9, f^{Ei,) and L^, and is independent of S, Y^, 
etc.. We plot the estimated values of rBs(^i/) determined 
from Eq. ([57]) in Fig. [TIT b). These indeed increase with 
L^. However, the estimated values of rBs{Lp) are al- 
ways larger than rx{Ly). This is because we have made 
many simplifications in deriving Eq. (|57p . In particu- 
lar, we have assumed that the "magnetic moments" of 
the NFIS's in the opposite directions are Sm?' /2{E^^) 
and 5m? /2{E^^) . This is a very crude approximation. 
According to Ref. [l3], the flavor conversion of neutri- 
nos and antineutrinos in the bi-polar mode is suppressed 
very little by the matter potential in the scenario with 
the inverted mass hierarchy. This is contrary to a con- 
temporary false belief that a large matter potential al- 
ways strongly suppresses neutrino flavor transformation. 
Therefore, rx{L„) should roughly trace the actual values 
of rBs{Li,), 

rx{L,) ~ rBs(i.). (58) 

We conclude that, like rBs(ii/), rxiL,^) with the inverted 
mass hierarchy has little dependence on S" or Yg, and in- 
creases monotonically with L^,. Because the single-angle 
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FIG. 12: (Color online) The difference between two almost 
identical systems grows exponentially in the region where col- 
lective neutrino flavor transformation changes from the syn- 
chronized mode to the bi-polar mode in the inverted neutrino 
mass hierarchy case. The solid line shows the exponential 
growth of |((5<;,/^)| (difference between the energy-averaged 
Ve_ NFIS's of the two systems) as a function of radius in 
the transition region. The dot-dashed line is a linear fit to 
In |(<5<;^5(r))|. This line has slope ~ 2.5/km. 



approximation uses the smallest value of rfjf among all 
trajectories, it gives lower bounds on the actual values of 
r^{L^). This is clear in Fig.[ll];b). 

We note that r-Qa,{L^) is the same for both the nor- 
mal and inverted mass hierarchies. Using this informa- 
tion, we can estimate whether significant neutrino flavor 
transformation for the normal mass hierarchy case be- 
gins in the bi-polar mode or not. Comparing panel (a) 
with panel (b) in Fig. [TTl we note that for the normal 
mass hierarchy, neutrinos and antineutrinos start flavor 
transformation in the synchronization mode for S — 250. 
They begin flavor transformation through the bi-polar 
mode for S — 140 when is less than a few times 10^^ 
erg/s. 

For the inverted neutrino mass hierarchy, our simula- 
tions show that chaotic behavior in neutrino flavor trans- 
formation can occur in the narrow region of radius where 
the collective behavior transitions from the synchronized 
to the bi-polar mode. To study this, we manually added 
random perturbations of order 10~^^ in ^j^^ [E) at radius 
r = 50 km (the region just below the synchronized-to-bi- 
polar transition for the particular case with S = 140 and 
= 10^^ erg/s). We follow the evolution of 5<;^{E), the 
difference between the NFIS's of neutrino v with energy 
E in the perturbed and unperturbed cases. We find that 
|(5i;y(i?)| for neutrinos of all species and energies grows 
with the same exponential factor in the transition region 



\5<;^{E)\ ocexp (2.5^ 



(59) 



In Fig. [m we plot ln|((5^y^)| as a function of radius r. 
Note that the difference T(^'>i'e)l between the energy- 



averaged Ve_ NFIS's of the perturbed and unperturbed 
cases is ~ 5 X 10* times larger at the radius where the 
system is fully in the bi-polar mode than in the region be- 
fore the synchronizcd-to-bi-polar transition. This chaotic 
behavior obviously causes difficulty in accurately simu- 
lating neutrino flavor transformation. However, we have 
performed several computations with different numerical 
schemes, all of which show qualitatively similar results. 
Therefore our analysis and conclusions are not affected 
by this behavior. At this point we do not know whether 
this behavior reflects true chaos or the appearance of a 
critical point in the neutrino/ antineutrino system. 

It is appropriate to comment on the validity of the 
single-angle approximation at this point. The traditional 
single-angle approximation picks the radial trajectory as 
the representative trajectory, which turns out to have the 
smallest n°*. As a result, the neutrino background effect 
tends to be underestimated. A slightly better approx- 
imation would be to average over all trajectories. This 
"averaged" case would have changed the geometric factor 
D{r/R^) in Eq. ^ to 



D{r/R.) = - 
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Even with this improvement, we would not expect to be 
able to simulate the complicated entanglement of neu- 
trino flavor transformation among different trajectories. 
The single-angle approximation is only accurate when 
neutrinos on different trajectories have the same flavor 
transformation histories. This does not seem to be the 
case for the bi-polar collective transformation. As a re- 
sult, it will be necessary to use multiangle simulations to 
accurately gauge, e.g., the effect of neutrino flavor trans- 
formation on Ye. Nevertheless, as we have demonstrated, 
the numerical simulations using the single-angle approx- 
imation are very useful as a means of exploring the ba- 
sic physics of neutrino flavor transformation in the hot 
bubble. These simple models do provide simple checks 
on more complex and computationally-intensive simula- 
tions. 



V. CONCLUSIONS 

We have carried out large-scale multiangle simulations 
of neutrino flavor transformation in the hot bubble em- 
ploying the atmospheric neutrino mass-squared differ- 
ence, \5m?\ ~ 3 X 10^"^ cV^, and effective 2x2 vacuum 
mixing angle 6 — 0.1. The numerical results we have 
presented support previous conjecture on the existence of 
collective neutrino flavor transformation of the bi-polar 
type in the supernova environment [iTj . Our simula- 
tions also show that both neutrinos and antineutrinos 
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can simultaneously undergo significant flavor conversion, 
largely driven by flavor ofl-diagonal potentials, at values 
of radius much smaller than those expected from ordinary 
MSW. This is along the lines of what was predicted in 
Ref. 16]. We have found that this flavor transformation 
occurs in both the normal and inverted neutrino mass 
hierarchy scenarios. 

For the normal mass hierarchy case, the full syn- 
chronization limit gives a lower bound on the radius 
where large-scale neutrino flavor transformation begins. 
(Ref. [T^ was the first to point out that, since i?sync is 
much smaller than average neutrino energies, synchro- 
nized flavor transformation modes operate closer to the 
neutrino sphere than those driven by the MSW matter 
potential.) Although an analytical analysis of neutrino 
flavor transformation in the bi-polar mode has yet to 
be done, our numerical simulations suggest that, for the 
normal mass hierarchy, the onset of bi-polar type flavor 
transformation always occurs at values of radius larger 
than those required in the full synchronization case. Our 
simulations also support the prediction of large-scale neu- 
trino flavor transformation in the inverted mass hierarchy 
scenario [13] • (Large-scale neutrino/antineutrino flavor 
transformation with small mixing angles in the inverted 
mass hierarchy case previously was seen in the early uni- 
verse context |32| and also in the supernova context [12] ■) 
We have found that this may occur at values of radius 
even smaller than those seen in the full synchronization 
mode in the normal mass hierarchy scenario. We have 
found that single-angle simulations can be used to give a 
lower bound on the radius where large-scale neutrino fla- 
vor transformation occurs in the inverted mass hierarchy 
scenario. 

Our "multiangle" calculations are the first to include 
self-consistent flavor evolution history entanglement on 
intersecting neutrino world lines. Although we flnd that 
"single-angle" simulations in some cases can give the cor- 
rect qualitative features of large-scale neutrino and/or an- 
tineutrino conversion in the late-time, hot bubble region, 
our simulations clearly show that a quantitatively cor- 
rect treatment must include coupled flavor development 
on different neutrino trajectories. Furthermore, since the 
location where large-scale neutrino and/or antineutrino 
flavor transformation begins in the supernova envelope 
can be a crucial issue for supernova shock reheating [ij, 
r-process nucleosynthesis ^ ^3 JTlJll [11 [ll , and the 
supernova neutrino signal 29, , is essential that 

simulations be quantitatively as accurate as possible. 

The simulations we have presented focus on the late- 
time supernova environment, i.e., the regime after the 
shock has been somehow re-energized. This epoch is a 
leading candidate for the site of the production of some 
or all of the r-process elements and will be a major fo- 
cus of future neutrino detectors/observatories should we 
be lucky enough to catch a galactic core collapse event. 
Though our simulations show that large-scale neutrino 
and antineutrino flavor conversion can take place during 
this epoch for the expected conditions of neutrino flux 



and entropy, we must go further than we have in this pa- 
per to produce quantitative predictions. There are three 
principal reasons for this: (1) we do not as yet know 
the matter density distribution above the proto-neutron 
star to sufficient accuracy at any epoch; (2) we do not 
know precisely the neutrino and antineutrino energy dis- 
tributions and fluxes which are emergent from the proto- 
neutron star; and (3) the matter composition (i.e.. Ye) 
can be affected by any changes in the neutrino and an- 
tineutrino spectra engendered by flavor transformation 
and we have not put this feedback in the calculations 
presented here. 

On point (2), recent work on supernova models at 
tpB < 1 s suggests that additional channels for neu- 
trino scattering may weaken or dilute the effects of the 
charged current opacities [H, [13]. This would tend to 
make the Vg, v^, and Vf^v^VrVT energy spectra more sim- 
ilar. Of course, if the neutrino energy spectra and fluxes 
are identical for all flavors, interconversion of these will 
have no astrophysical effect. We note, however, that re- 
liable neutrino transport calculations at the late time we 
considered here do not exist, and the core's composition 
and neutron excess is expected to change considerably 
between tpB ~ 1 s and tpB — 10 s. Clearly, this issue 
is critical for gauging the astrophysical effect of neutrino 
flavor mixing. 

It is well known that density fluctuations on short 
length scales and other inhomogeneities can modify co- 
herent neutrino flavor evolution through MSW reso- 
nances (35l . [36| . How these fluctuation-induced modi- 
fications could manifest themselves in quantum flavor 
history entanglement on intersecting neutrino trajecto- 
ries is not known. This issue may be closely related to 
the problem of calculating neutrino transport and pre- 
dicting the emergent neutrino energy spectra in general 
[37, 38, 39, 40, 41] and to the inclusio n of neutrino fla- 
vor mixing in the core in particular (4^ . [43j . Though 
our simulations are spherically symmetric, they do show 
that the density and Yg profiles near the proto-neutron 
star surface are important for obtaining the correct fla- 
vor evolution of the neutrino and antineutrino fields, even 
well above the proto-neutron star. 

These uncertainties aside, our calculations indicate 
that large modifications of the emergent neutrino and an- 
tineutrino energy spectra are likely to occur over most of 
the range of expected thermodynamic and neutrino emis- 
sion parameters of relevance in the late-time supernova 
environment. Furthermore, we have found that these 
modifications could set in sufficiently deep in the super- 
nova envelope to affect Ye tE] and r-process nucleosyn- 
thesis [13, [4^ through neutrino interactions. However, 
we have not included charged-current weak interaction 
]Eq. (|¥T|l ] feedback in the calculations presented here. 

We have fixed Ye and gs in this work, essentially to 
simplify the computations. In future simulations we will 
remove these constraints and allow Yg, and gs to be cal- 
culated consistently with feedback from neutrino capture 
reactions. However, we expect that the collective neu- 
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trino flavor transformation illustrated here will not be 
changed qualitatively with changing Y,, and gs- The bi- 
polar neutrino flavor transformation seen in our simula- 
tions is largely independent of the values of Ye and gs- 
For example, we have shown that i>e of energy smaller 
(larger) than a critical energy Eq could convert to other 
flavors if Sm^ > {6m^ < 0). This critical energy Eq 
asymptotically approaches a limit if L^, is large enough, 
or equivalently, the electron density profile is sufliciently 
condensed toward the proto-neutron star. The asymp- 
totic limit of Ec depends only on the neutrino mixing 
parameters and the initial energy spectra for neutrinos 
and antineutrinos. 

Because the proto-neutron star is neutron-rich, the ini- 
tial energy spectrum may be softer than those for neu- 
trinos in other flavors. Our simulations suggest that I'e 
and neutrinos in other flavors may swap the low-energy 
{El, < Ec) or high-energy {E„ > Eq) parts of their 
spectra depending on the sign of ^771^3. Note that this 
stepwise swapping is independent of the details of the 
neutrino energy spectra. With other effects correctly ac- 
counted for and a good signal from a galactic supernova, 
this phenomenon may offer a unique probe of the neu- 
trino mass hierarchy problem. 

We have employed 2x2 neutrino flavor mixing in our 
simulations. It is possible to extend our codes to im- 
plement the neutrino mixing of all three active flavors. 
However, we expect that neutrino flavor transformation 
in the hot bubble region will not change much on in- 
clusion of a third neutrino flavor. For one thing, 
and Vr are almost equally mixed in the hot bubble be- 
cause they experience the same weak interactions and 
023 — 7r/4. For another, the two neutrino mass-squared 
differences, Sml^j^ and Sm^, are separated by over an 
order of magnitude. Taking Sm? = Srn^ ~ 8 x 10~^ eV^ 
and 6 = 612 ~ 0.6, we estimate that the onset radius 
of large-scale neutrino flavor transformation in the full 
synchronization limit is 7'x/MSw(-E'sync) — 227 km for an 
entropy per baryon S = 140. This location is almost 
outside the range of the simulation results presented in 
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In summary, though many aspects of our calculations 
are reasonable approximations at best {e.g., 2x2 mix- 
ing, assumptions of spherical symmetry, an inflnitely thin 
neutrino sphere, neutrino/ antineutrino energy spectra of 
the Fermi-Dirac type, etc.), our computations do mark 
an important advance in that they self-consistently treat 
coupled neutrino flavor evolution on different trajecto- 
ries. We cannot claim generality for our conclusions. 
However, our assumptions are reasonable, and our re- 
sults are robust, and so there is nothing to suggest that 
our results represent an isolated case either. Not only 
do our results show that a proper treatment of coupled 
neutrino trajectories is important, but they also indicate 
that the measured neutrino mass-squared difference val- 
ues and mixing angles likely imply large-scale flavor con- 
version of neutrinos and antineutrinos in astrophysically 
important regions in the post-explosion supernova envi- 
ronment. 
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